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Three-dimensional orientational order in systems whose ground states possess non-zero, chiral 
gradients typically exhibits line-like structures or defects: A lines in cholesterics or Skyrmion tubes 
in ferromagnets for example. Here we show that such lines can be identified as a set of natural 
geometric singularities in a unit vector field, the generalisation of the umbilic points of a surface. 
We characterise these lines in terms of the natural vector bundles that the order defines and show 
that they give a way to localise and identify Skyrmion distortions in chiral materials - in particular 
that they supply a natural representative of the Poincare dual of the cocycle describing the topology. 
Their global structure leads to the definition of a self-linking number and helicity integral which 
relates the linking of umbilic lines to the Hopf invariant of the texture. 


Three-dimensional systems described by orientational 
order, and more generally vector fields, are often seen 
to contain line-like geometrical features. For example, 
in cholesteric liquid crystals as A lines im or disloca¬ 
tions, and as double twist cylinders in blue phases Ed, 
in superfluids they arise as the cores of vortices mm, 
and they can be identified in the Skyrmion [12] textures 
of chiral ferromagnets mm and Bose-Einstein conden¬ 
sates [T8K23] . These structures are being manipulated 
with increasing precision in experimental systems, pro¬ 
viding a variety of topologically non-trivial textures in 
ordered media with ample potential for the development 
of structured metamaterials [241126] , fluidic photonics m 
and ultralow current spintronics [281130] . Skyrmions in 
chiral magnets are often described as magnetic whirls 
that encode a topological winding. The same topological 
profiles can be imprinted on spin configurations of Bose- 
Einstein condensates, superfluids and electronic states. 
Skyrmions also arise in liquid crystals [DISIHSS], where 
the nature of nematic order allows for a greater variety 
of three-dimensional textures including knots, a focus of 
activity across several disciplines [36d4T] , which can then 
be entangled by Skyrmions [42] . All of these display a 
blend of topological features, giving robustness, with a 
geometrical character that comes from energy minimisa¬ 
tion. In this paper we show that these objects can be 
described as a set of natural geometric singularities in 
a vector field, which we call umbilic lines, whose topo¬ 
logical significance comes through their representation of 
the Poincare dual to the Euler class of a naturally defined 
vector bundle. 

The umbilics we describe are generic traits, exhib¬ 
ited by any vector or line field, and closely analogous 
to several other structural degeneracies, for instance the 
lines of circular polarisation (C lines) in a generic elec¬ 
tromagnetic field [43H46] and the umbilic points of sur¬ 
faces mm and random fields [49U56] . from which the 
terminology derives. Indeed, degeneracies in geometri¬ 
cal physical properties are commonplace; their relevance 
and significance for the description of global or struc¬ 
tural aspects of physical behaviour, and the identifica¬ 


tion of generic traits, has been cemented since Berry’s 
influential paper on the geometrical phase in quantum 
mechanics mm- These degeneracies are not always 
random; for instance, it has been known for a long time 
that the clear blue sky is polarised and that there are 
points of degeneracy (or singularity) in that polarisation, 
the Babinet and Brewster points either side of the sun, 
and the Arago and second Brewster points either side of 
the anti-sun. Away from these points the polarisation 
is elliptical and the degeneracy is to a state of circular 
polarisation, about which the principal axes of the po¬ 
larisation ellipse rotate by 7r, corresponding to a +1/2 
index. That there are four such points is a reflection of 
the fact that the sky has the topology of a sphere and an 
illustration of the Hopf index theorem 0HEDI- The same 
Berry phases control aspects of the response in magnetic 
Skyrmions [6TH63] . 

We give general local and global descriptions of ori¬ 
entational order in three dimensions, identifying umbilic 
lines as a key natural feature characterising the ordered 
material. The flux of umbilics through any surface is 
identified with the Skyrmion charge through an applica¬ 
tion of the Gauss-Bonnet-Chern theorem. The umbilics 
convey the topology in the orientational order since they 
are Poincare dual to twice the Euler class of the two- 
plane bundle of directions that are everywhere orthogo¬ 
nal to the local director. The extended nature of umbil¬ 
ics entails additional global properties that we describe 
in terms of self-linking numbers. An analogue of helicity 
connects to the Hopf invariant in a manner parallel to 
that for fluid flows [641166] . Of additional note is a con¬ 
nection of our local description with recent work of Efrati 
and Irvine m on chiral shapes and structures. In their 
work they introduced a chirality pseudotensor to char¬ 
acterise the direction of twisting and sense of chirality 
of various materials. We demonstrate a general connec¬ 
tion with the shape operator for a surface or vector field, 
from which their chirality pseudotensor can be derived, 
and show that natural benefits are gained by considering 
both operators. Analogous geometrical structures to the 
umbilics we describe have been discussed by Thorndike et 
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al. [68], Copar et al. [69], and Beller et a 1. [70], although 
not with the topological aspect that we give here. 

Umbilic lines acquire extra significance in chiral ma¬ 
terials where the ground state has non-zero gradients. 
In these materials the vortex configurations surrounding 
umbilics correspond to local minima of the free energy. 
For concreteness, then, we consider ordered materials de¬ 
scribed by a unit magnitude vector field n, called the 
director, and whose free energy has a form similar to 


F 



+ 2g 0 n • V x n 


d 3 x, 


(i) 


where K is an elastic constant and qo is the strength of 
the chirality, or magnitude of spin-orbit coupling. This is 
the Frank free energy m for a bulk cholesteric under the 
one elastic constant approximation and also the Hamil¬ 
tonian for elastic distortions in a chiral ferromagnet with 
a Dzyaloshinskii-Moriya interaction dies]. 


LOCAL GEOMETRY OF ORIENTATIONAL 
ORDER 

The local geometry of the orientational order at any 
point in the material can be given the following generic 
characterisation. Since n is a unit vector its gradients 
are always perpendicular to itself. They separate natu¬ 
rally into gradients along the direction of the local ori¬ 
entation, known as bend distortions, and gradients along 
orthogonal directions, which correspond to splay, twist, 
and saddle-splay deformations. These latter represent a 
linear transformation on the two-dimensional plane of di¬ 
rections perpendicular to the local orientation, which is 
how we shall think of them. If di, d 2 are unit vectors or¬ 
thogonal to n, such that {n, di, d 2 } form a right-handed 
set, then the orthogonal gradients can be expressed in 
this local basis as a 2 x 2 real matrix with the decompo¬ 
sition 


V_Ln 


"1 0" 

n • V x n 

"0 -1 

0 1 

+ 2 

1 0 


Ai A 2 
A 2 —Ai 


( 2 ) 


The first component is an identity transformation on the 
orthogonal directions, while the second is a skew trans¬ 
formation that provides the orthogonal directions with a 
complex structure. The coefficients of these two terms 
are known as the mean curvature and the mean tor¬ 
sion [ZH [ZE], respectively, and correspond to scalar and 
pseudoscalar parts of the transformation. The last part 
is a deviatoric, or spin 2, component, A, that conveys 
the local ‘shear’. Its eigenvectors can be used to define 
principal directions of the curvature of the local orienta¬ 
tion. In the case where n is the normal to a surface, the 
transformation © is the derivative of the Gauss map, 


or shape operator of the surface, and the principal direc¬ 
tions alluded to are the directions of principal curvature. 

The geometric character of the orientational order is 
conveyed by this generic local description. For instance 
the twist, n • V x n, encodes the chirality of the order 
with surfaces of zero twist separating regions of different 
handedness. Zeros of the deviatoric part of V^n occur 
along lines in three dimensions - they require both Ai 
and A 2 to vanish and hence are codimension two - and 
correspond to degeneracies of the principal curvatures, or 
singularities of the principal curvature directions. They 
are the umbilics of the orientational order. At such points 
the transverse gradients are isotropic meaning that um¬ 
bilics are naturally associated with the centres of vortex- 
like structures. A basis-independent expression can be 
given for these zeros purely in terms of the director; they 
correspond to the vanishing of the non-negative quantity 


|A | 2 




n • V x n 


n(V • n) - (n • V)n , 


(3) 


this being the 2-norm of the matrix A. In the vicinity of 
an umbilic we may write 


A = |A| 


cos 6 sin 0 
sin 6 — cos 0 ’ 


(4) 


and around the umbilic the angle 0 winds by an inte¬ 
ger multiple of 27r. An umbilic is said to be generic if 
it has a winding number of =bl. Umbilics with higher 
winding do occur, and indeed those with winding +2 
are often observed in vortices, double twist cylinders and 
Skyrmions, with the degeneracy coming from an imposed 
axisymmetry. This can persist in certain experimental 
settings, either where there is axisymmetry or hexagonal 
symmetry, however, from a structural point of view they 
can be seen as a composite of winding ±1 umbilics, into 
which they break apart under perturbation. The generic 
umbilic points of a surface have been classified as stars, 
lemons, or monstars [48], according to the local profile 
that the principal curvature directions have in the vicin¬ 
ity of the umbilic point. The same local structure arises 
in the rotation of the eigenvectors of A around the um¬ 
bilic lines of a generic vector field in three dimensions - 
it is not hard to see that the eigenvectors of A rotate by 
only half as much as A itself, i.e. by a half-integer mul¬ 
tiple of 27r, this being the familiar phenomenon of a real 
Berry phase [571158] - so that the local profile of an um¬ 
bilic line may be classified in the same way as the umbilic 
points of a surface or the C lines of generic electromag¬ 
netic fields [43]. As this description is classical |48] [68] 
we choose to focus on complementary aspects, making 
more use of the language of vector bundles m • 

Illustrations can be given for a number of simple field 
conformations corresponding to textures observed in su¬ 
perfluid vortices, magnetic Skyrmions and the A lines 


















3 


emeses cacao «=><£><?<£> <9 <j <y & <y o o o 
QSQQQOO 0^(?o 0 O O O O €> €> €> 
CCCCWO& q 0 0 0&&€>&&& 
oooe3c5tf5^ o 0 0 0 ##<?€>&&& 
oarac=c=ie3<JJ> 0 0 0 0 0 0€>€?t2'e3'S3 

ooooeaeac? oODB6<*&<* c * c > ct 

OOOOOO0B 00g@o 00000 

0000®*3*)$ftft% 

»»■>«)«,%,% o^H 

cananocs^ o ^ ^ ^ ^> <S<S< 55 C 5 C 3 

C3C3C3«=5C3C3C3 %) 0 O^^^QQOCSO 

#<*«*O00O<S>%fcO6C>CCOOC) 

raraesesesCaSa C3 <S<S<S 


C 


W+OOOVOM** 



ftftft 
ft ft ft 

oo 

ft ft ft 
ft ft ft 
ft D ft 
D D D 
0 D 0 
ODD 
D D D 

o a o 
a o o 

ODD 

ODD 

DDD 


ftftftftftftftft 
ftftftftftftftft 
ftftftftftftftft 
ftftftft^ftftft 
ftft^^^ft 
ftft<S<S^<Sft 
ft ft<s^>^^><sft 
D []*-,====«] 
D0#&&&#0 
0 0#<?<?<?# 0 
0 0 00000 0 
00000000 
00000000 
00000000 
00000000 
00000000 


ftftftftft 
ftftftftO 
ftftftftQ 
ftftftftO 
ftftftftft 
ftftftftO 
ft ft 0 ft D 
ft ft ft D D 
D D D D □ 
DO DDD 
ODD 0 0 
00 D 00 
0 0 D D 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 D D 




FIG. 1: Umbilic lines, (a) A + A - dislocation in a cholesteric. 
It contains two generic umbilics. This combination of two 
generic umbilics can be identified as half a Skyrmion, or a 
meron. (b) Eigenvector field of II for the dislocation. The 
two umbilics have opposite windings of the eigenvectors, but 
are counted with the same sign as n rotates by 7r when pass¬ 
ing from one to the other, (c) Umbilic lines can be oriented 
in three dimensions, (d) & (e) Typical centres of vortices, 
based on curvature and torsion distortions respectively, that 
correspond to winding number 2 non-generic umbilics. 


of cholesteric liquid crystals, shown in Fig. [l] Two- 
dimensional sections of the full textures are shown on sur¬ 
faces intersecting the umbilics transversally. The eigen¬ 
vectors of II shown in (b) identify the profiles of the um¬ 
bilics in (a). Away from the umbilics this eigendirection 
coincides with the pitch axis of the helical texture. 

We explain later that the windings just described are 
equivalent to the indices of the umbilics up to a sign, 
which arises because the orientation of n reverses from 
one umbilic to another. (The discrepancy is therefore 
absent if n is the normal to a surface, since then there 
is no reversal of its orientation between umbilics, and 
the index is simply the local winding of the eigenvectors 
of A.) The sum of the indices of the umbilics yields a 
topological invariant, equal to four times the Skyrmion 
number of the texture on the two-dimensional surface. 
This calculation of the Skyrmion number is dual to the 
usual method of integrating a curvature form H2ES1 

q = — [ e a i) C n a din b djn c dx l A dx\ (5) 

which computes the degree of the map n : E —)> S' 2 . 
The duality can be viewed as entirely analogous to that 
in the classical differential geometry of surfaces between 
the Gauss-Bonnet theorem - integral of the Gaussian cur¬ 


vature - and the calculation of the Euler characteristic 
using umbilic points of the surface. 

Two additional aspects of the local description are 
worth bringing out before giving a global treatment. 
First, the complex structure J = [^“o 1 ] acts on the 
shape operator © by simple composition to produce an 
equivalent linear transformation 


T __ n • V x n 

"1 

o' 

Vn 

'0 -1 

x = JoVin = 

_0 

1 

+ 2 

1 0 _ 


—A 2 Ai 
Ai A 2 


(6) 


This transformation, although equivalent to the shape 
operator, is independent from it in the sense that they 
are orthogonal with respect to the natural inner prod¬ 
uct between matrices, Tr(x T V^n) = 0, a property that 
turns out to be of some importance in the analysis of the 
topology of umbilics. x is the chirality pseudotensor in¬ 
troduced recently by Efrati and Irvine [671 i n their study 
of chiral materials; our description of it provides a com¬ 
plementary perspective to theirs. The spin 2 component 
II = J o A carries equivalent information to A. Its eigen¬ 
vectors can be thought of as defining directions of prin¬ 
cipal torsion (or twist) and are simply rotated relative to 
those of A by 7r/4 about the director. Nonetheless, for 
materials with chiral interactions or structure, it is these 
directions, rather than the eigenvectors of A, that are of¬ 
ten more visibly recognisable, as is elegantly illustrated 
in the experiments of Armon et a 1. ca on the opening 
of chiral seed pods, and those of Efrati and Irvine [67] 
on stretched elastic sheets. For this reason, we plot the 
eigenvectors of FI rather than those of A in all figures in 
this paper. 

Second, if the umbilic forms a closed loop then the 
winding of the angle 6 along the contour length of the 
umbilic (a longitude) conveys a second integer invariant. 
Of course, this latter depends on both the choice of lon¬ 
gitude and the homotopy class of the local trivialisation 
{di,d 2 }. We describe this in detail later. 


GLOBAL DEFINITION OF UMBILICS AND THE 
TOPOLOGY OF VECTOR FIELDS 

The umbilic lines we have identified convey topological 
information about the vector field n, specifically the Eu¬ 
ler class of the orthogonal 2-plane bundle, £, through a 
combination of the Gauss-Bonnet-Chern theorem [H6H83] 
and Poincare duality. The curvature of a vector bundle 
represents a characteristic cohomology class that can be 
integrated over homology cycles to produce topological 
invariants (Euler numbers). By Poincare duality this can 
be viewed as a homology class and the umbilics furnish 
a representative of this Poincare dual. We shall develop 
this general topology of orientational order in terms of a 
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FIG. 2: The splitting TR 3 = L n ® £ induced by orientational 
order. Left: A simple helical director profile, corresponding to 
the ground state of a cholesteric liquid crystal or helimagnet. 
Right: An idealised, axisymmetric Skyrmion profile, the cen¬ 
tral portion of which is a standard idealised model of double 
twist cylinders in blue phases. 


natural connection on the vector bundle that the spin 2 
transformation A takes values in. 

As before, let n be a unit vector field in M 3 ; then it 
defines a splitting of the tangent bundle 

TR 3 “ L„ © (7) 


into a real line bundle, L n , of vectors parallel to n and a 
rank 2 real vector bundle, £, of vectors orthogonal to n. 
Examples of this splitting are shown in Fig. [2] The gra¬ 
dients of n take values in TM 3 * 0 £ = (L* 0 £) ® (£* 0 £), 
with the first part containing the bend deformations and 
the second the perpendicular gradients given locally in 
By complete analogy with the classical differential 
geometry of surfaces we call these perpendicular gradi¬ 
ents the shape operator of the orientational order and 
think of it as a linear transformation on the orthogonal 
2 -planes. 

The shape operator decomposes in the following way. 
Local rotations about the director field correspond to a 
local symmetry and endow the bundle £ with a natu¬ 
ral SO( 2 ) action, under which the tensor product £* 0 £ 
splits into a direct sum X ® J ® E, where X and J are 
trivial line bundles and E is a rank 2 vector bundle. Cor¬ 
respondingly, the shape operator has the decomposition 


V_Ln 


Vn 

2 


k + 


n • V x n 

2 


J + A, 


(8) 


given previously in (J2|. Here, is the identity transfor¬ 
mation and J — nx is a complex structure on £, now 
defined globally, while A has the global definition 


A(v) = i (v • V)n + n x (n x v • V)n 


(9) 


for any v in £. The umbilics are identified with the zeros 
of A; they encode the topology of the orientational order. 

Without repeating too much of the local description 
given previously, we record simply that composition with 
the complex structure again yields the chirality pseu¬ 
dotensor m, X = J o Vj_n, and that the spin 2 com¬ 
ponent of this, n = J o A, can be expressed globally 


as 


1 r 


n (v) = 2 


n x (v • V)n — (n x v) • Vn 


( 10 ) 


for any v in £, and is also a section of E. We reiterate 
that although n and A may be obtained one from the 
other by composition with the complex structure, they 
are orthogonal with respect to the natural inner product 
(n, A) = ^Tr(n T A). At the same time, they share the 
same magnitude and the same zero set, the umbilics U. 
Thus, away from umbilic lines, n and A represent two 
non-zero orthogonal sections of E, and so they provide 
a basis for E\ R s_ u and also define a natural connection 
1 -form 


„ _ <n,VA> 
(A, A) • 


( 11 ) 


This 1-form conveys the fundamental topology of the ori¬ 
entational order. 

We give in ( 11 ) a global definition of A that emphasises 
its naturality and makes clear how it arises. In a local 
chart it can be represented in terms of the eigenvectors, 
p + and p - , of A and has components A{ = 2p~dip+, 
an expression perhaps closer in form to the usual way of 
writing Berry connections. However, one must keep in 
mind that this representation is only local in the present 
context; otherwise it only serves to obfuscate the natural 
structure. The exterior derivative of A, LIe = dA, is the 
(Berry) curvature of the vector bundle E. It is equal to 
twice the curvature of £ on account of E being a spin 2 
subbundle of £* 0 £. (This is also evident from the local 
expression in terms of the eigenvectors p ± .) 

A fundamental property of an umbilic is given by the 
integral of the connection 1 -form around any simple 
closed loop 7 encircling the umbilic. This is the holonomy 
of the connection, or Berry phase. Let E be an oriented 
surface with boundary 7 , as in Fig. HI If s is generic 
then the umbilics will intersect it transversally in some 
number of distinct points. Let {Di} be disjoint small 
discs in E about each of these intersection points and 
denote their boundaries by {Ci}. Then |A| is nowhere 
zero on E — {Di} and d(E — {Di}) = 7 U {—Ci} so that 
by Stokes’ theorem 


/[ A = [ dA = / 

j 7 7 JCi Jt. 


LI F 




( 12 ) 


To extend the analysis over the umbilics, we can intro¬ 
duce a local trivialisation of L n ® £ on each of the Di, 
which we denote by the orthonormal basis {n, di,d 2 }. 
Then ei = di 0 di — d 2 0 d 2 and e 2 = di 0 d 2 + d 2 0 di 
form a local basis for E. In such a basis A = | A| cos 6 e 1 + 

| A | sin#e 2 and A = d6 + 2uj where 0 di = (d 2 ) a <9^(di) a is 
the connection 1 -form on £ for this basis (the factor of 
two accounts for E being a spin 2 bundle). With this, 
one finds at once the Gauss-Bonnet-Chern theorem in 
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the form 



= + index s U u 


(13) 


which defines the index of an umbilic, an integer associ¬ 
ated to its intersection with a surface. The sum of these 
over all of the umbilics is an Euler number of the bundle 
E. Note that, as alluded to earlier, the index of an um¬ 
bilic is not simply the winding number of the angle 6 that 
characterises the local profile as that is measured relative 
to a local basis for £, while the index uses the orientation 
on the surface E; there is, therefore, an additional sign 
according to whether dO is cooriented with Ci or not. 

It is instructive to contrast this situation, where the 
type of profile does not carry topological information, 
with that of umbilical points of surfaces where the type 
of profile does carry topological information about the 
sign of the winding. This occurs because in the case of 
surfaces n is fixed to be the normal to E, and so dO is 
always cooriented with Ci. 

One can understand this relation between the profile 
and index in terms of oriented umbilics. The 1-form A 
circulates around umbilic lines, orienting them such that 
the circulation is counter-clockwise. If one orients um¬ 
bilics in this way, the expression for the index in (13) is 
rewritten as 


indexs Ui = sJnt(L^,E), (14) 


where Int(L^, E) is the signed number of intersections of 
Ui with E and Si E N is the absolute strength of the um¬ 
bilic Ui. Examples illustrating this are shown in Figs, [l] 
and [4] we describe only the latter (the former is essen¬ 
tially identical). Fig. [ 4 ] shows a Skyrmion lattice, a well- 
known structure seen in a variety of systems m lua- 
HZl ESI ESI I84H86] . The umbilics form a lattice, with 
8 = 2 umbilics each surrounded by six 8 = 1 umbilics. 
The eigenvectors of II show the central umbilic having a 
+1 profile, and the six outer umbilics a —1/2 profile. All 
the umbilics have the same orientation induced by A, and 


so their intersection numbers, and thus indices in (13), 


are all positive. For the +1 profile, —1/2 profile and the 
orientations to all be consistent, it is necessary that n 
rotates by 7 r when passing from the central umbilic to an 
outer one. 

The curvature of the bundle £ can be expressed in 
terms of the director field as 


= —-eaijc^dir^djjfdx 1 A dxE 


(15) 


This is because each 2-plane of the bundle £ is explic¬ 
itly embedded in M 3 (see Fig. [ 2 ]) and the expression is 
just the extrinsic definition of curvature in terms of the 


A\/xB 



FIG. 3: Anatomy of an umbilic line. Top. The vector V X 
B is tangent to generic umbilic lines, giving them a natural 
orientation. This orientation is compatible with that induced 
by A. Middle. The vector field n and associated orthogonal 
plane £. Note that in general the tangent vector of the umbilic 
does not align with n. £ is oriented by n and so gives a 
natural cycle along which to measure the winding. For a 
generic umbilic if n • V x B > 0 then the eigenvectors of A 
have a +1/2 profile and if n • V x B < 0 the eigenvectors have 
a — 1/2 profile. Bottom. A measuring cycle 7 bounding a 
surface E. Whether 7 measures the winding of the umbilic as 
positive or negative depends on whether the umbilic - oriented 
by V x B - intersects E positively or negatively. 


(generalised) Gauss map. Equally, it follows from an ex¬ 
plicit calculation in local coordinates for any particular 
texture, say a radial hedgehog. A third way to establish 
it is through a calculation in the style of the Mermin-Ho 
relation [ 8 ], a good account of which has been given by 
Kamien m • The curvature of E is equal to twice this, 
Qe = since E is a spin 2 subbundle of £* ®£. There¬ 
fore, when E is closed, or the boundary term vanishes, 
we have 

sJnt(£/i,E) = [ e a i) C n a din b djn c dx l A dx j = 4g, 

i 27r Je 

(16) 

or the sum of the indices of the umbilics piercing E is 
equal to four times the Skyrmion number. This estab¬ 
lishes a general relationship between Skyrmions and um¬ 
bilics. This is illustrated in Fig. [4] where each unit cell 
contains one 8 = 2 umbilic and six 8 = 1 umbilics, each of 
which are shared between three unit cells, giving a total 
count of 2 + 6 x 1/3 = 4, i.e. there is one Skyrmion in 
each unit cell. This approach can also be extended to the 
charcterisation of merons, or half-Skyrmions, which can 
be thought of as a single 8 = 2 umbilic, or two 8 = 1 um¬ 
bilics, as in the case of the cholesteric dislocation, shown 
in Fig. [l] 

This correspondence between Skyrmions and umbil¬ 
ics can be seen as a version of Poincare duality, giving 
a natural way to localise Skyrmions to a set of points 
or lines. The set of umbilics U in a three dimensional 
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FIG. 4: Skyrmion lattice. Each unit cell possesses a single 
central s — 2 umbilic, with the six outer s — 1 umbilics each 
shared between 3 unit cells giving a total count of 2+6 x 1/3 = 
4 umbilics = 1 Skyrmion per unit cell. Top Left : Unit vector 
field n in a Skyrmion lattice. Top Right : Contour plot of |A| 
showing the umbilic lines forming a lattice of tubes. Bottom 
Right : Eigenvector field of II, illustrating the +1 winding 
around the central umbilics and —1/2 winding around the 
outer umbilics. Note that both are counted with the same 
sign as the orientation of n changes from one to another. 
Bottom Left : Umbilics form lines in three dimensions. 



FIG. 5: Transient monopoles and umbilic lines observed in 
a simulated quench of the free energy 0. Left: Contours 
of constant value of |A|. The blue lines are umbilics, the 
beige spheres are monopoles (point defects). Each umbilic 
has winding n = 1. They are all oriented similarly, so that 
the sum of the umbilics entering/leavi ng th e negative/positive 
point defect is 4, in accordance with |l3| and corresponding 
to the Skyrmion created between the monopoles. Top Right: 
Cross-section showing n on the shaded plane. As a map from 
the plane to S' 2 , it has degree 1, indicating the presence of a 
Skyrmion. Bottom Right: Contour plot showing the magni¬ 
tude of log(| A|) on the shaded plane, the umbilics are clearly 
picked out. 


contains a Skyrmion, illustrating the well-known creation 
and annihilation of Skyrmions in terms of point defects, 
observed experimentally in m 


domain M gives a representative of a cycle in 
The Poincare dual of this is the cocycle e(E) E JT 2 (M), 
the Euler class, which depends only on the topology of 
the bundle E and hence n, and gives a global constraint 
on the total number of umbilic lines in the system in 
terms of the topology of n. As a spin-2 vector bundle, 
e(E) = 2e(£) = 4[n], where [n] E H 2 (M) is the cocy¬ 
cle represented by n [88]. In the absence of torsion in 
H 2 (M) [89], there is a canonical representative of the co¬ 
cycle 4[n] as a 2-form which is given through the Chern- 
Weil homomorphism by Pf(—^/2 tt), which leads to (13). 

An example of the Euler class as a global constraint 
on umbilics is given by the case of monopoles, illustrated 
in Fig. [5] These are point defects in the magnetisation 
that mediate changes in the number of Skyrmions El- 
In the local neighbourhood of a point defect the topol¬ 
ogy of n is described by an integer q E ^(S 2 ) = Z, the 
charge of the point defect. Through (Jl6| one finds that a 
sphere surrounding this point defect must be pierced by 
umbilics of total index 4 q. In Fig.[5]we show a simulation 
of two ±1 point defects, or monopoles, in a chiral ferro- 
magnet. There are four umbilics, each of generic type, 
ending on the point defects in a manner reminiscent of 
Dirac strings EH- On a slice in the midplane of the cell n 


EXAMPLES OF LINE FIELDS: BLUE PHASES 

So far we have confined our analysis to those systems 
where the orientational order is represented by a vec¬ 
tor. In liquid crystalline phases the orientational order is 
instead line-like. Perhaps the most noticeable feature 
of this is the presence of line defects, called disclina- 
tions, in liquid crystals [90, 91 , which make manifest the 
non-orientability and are absent in materials with vector 
order. However, they can arise even in vector ordered 
materials if the vector couples to an internal degree of 
freedom to give nematic-like symmetry; examples include 
thin films of 3 He [92] and spin 1 polar condensates [93 . 

As the director is non-singular along umbilics, the ef¬ 
fect of the non-orientability is primarily restricted to their 
global identification. Indeed, since umbilics can be iden¬ 
tified from the local structure of the director gradients, 
they have the same local description for line fields as they 
do for vector fields. One need only choose an orientation 
for the line field, which may always be done locally, i.e. 
in regions that do not contain disclination lines. There 
is even the same description for their global character 
in that the umbilics still furnish a representative of the 
Poincare dual to twice the Euler class of the orthogonal 
2-plane bundle £. Interestingly, for line fields with knot- 















7 


ted disclination lines the Euler class of the orthogonal 
2-plane bundle is usually torsion [42], which the usual 
Chern-Weil theory of curvature forms is insensitive to. 
(The curvature form represents only the integral part of 
the cohomology [76l 153] .) The integral in (16) thus does 
not provide any information in such cases. However, as 
the umbilic lines are Poincare dual to twice the Euler 
class of £ (and 4[n]) they continue to capture the topol¬ 
ogy of the orientational order, even when it is torsion, 
modulo elements of order 4. 

As discussed earlier, and shown in Fig. [lj one can iden¬ 
tify umbilics with A lines in cholesterics. In this way, 
(13) can be seen as a global constraint on the A lines in 
a cholesteric, as determined by the topology of n, with 
additional subtlety associated to elements in H 2 (M) of 
order two and four. We anticipate that the elements of 
order 4 are associated with r lines. For the elements of 
order 2 the orthogonal 2 -plane bundle £ admits a global 
section and we anticipate that these are associated with 
distinct “cholesteric ground states” m- 

Rather than present a detailed account of the anal¬ 
ysis of line fields, we focus on simple illustrations and 
a practical means for identifying umbilics globally, us¬ 
ing one particular physical system - the cholesteric blue 
phases. These remarkable materials are characterised by 
periodic orientational order corresponding to crystalline 
space groups, while remaining three-dimensional fluids. 
Two distinct cubic structures, BPI with space group 
0 8- (/4 i 32) and BPII with space group 0 2 (P4 2 32), are 
observed on varying temperature or the amount of chi¬ 
ral dopant 0131251, while several other structures can be 
induced by applied electric field effects l96ff98l or confine¬ 
ment [99] . and there is also an amorphous phase, known 
as BPIII |100| . Their structure arises because the lo¬ 
cal optimal configuration for the chiral free energy is one 
of double twist (illustrated in Fig. [l] and [2]) , but this is 
only locally favourable and the structure becomes ener¬ 
getically unstable if it extends too far. This leads to a 
frustrated arrangement where locally preferred regions 
of double twist form periodic arrays interwoven by a net¬ 
work of disclination lines urns. The blue phases have a 
long-standing association with Skyrmions in chiral mag¬ 
nets, with which they may be considered analogous. This 
analogy continues to stimulate considerable interchange 
of ideas between the two fields njEmonm 

It is usual to describe the order in liquid crystals, es¬ 
pecially in numerical simulations, using a tensor order 
parameter, called the Q-tensor, that captures the non- 
orientability of the director field m- The Q-tensor is 
proportional to the deviatoric part of the dielectric ten¬ 
sor or the second moment of the molecular orientational 
distribution function. The director field can be recov¬ 
ered from it as the eigenvector associated to the largest 
eigenvalue. In typical nematic materials, the Q-tensor is 
uniaxial and we have Q oc n(g)n— |l. Our first objective 
is to describe how umbilics may be identified from the 


Q-tensor. 

Orientational order continues to define a splitting of 
the tangent bundle TM 3 = L n ®£ whether it is orientable 
or not, except that the line bundle L n is non-trivial in the 
non-orientable case, with the consequence that the direc¬ 
tor n is not globally defined (as a vector). Rather, the 
projector n0n provides a section of L n <g> L n . Similarly, 
the shape operator has to be modified so as to be pre¬ 
sented as a section of the twisted bundle £* m- 

However, interestingly, the chirality pseudotensor, and 


its deviatoric part n ( 10 ), continue to be globally well- 


defined in the same form even when the director is non- 
orientable, so that n admits an immediate extension to a 
‘Q-tensor version’. We shall use this to identify umbilics 
in line fields. 

As £* (8) £ is a subbundle of TM 3 * ® TM 3 the linear 


transformation n, equation ( 10 ), can be embedded in a 


3x3 matrix and expressed, in a Cartesian basis, as 


n ij — t^ilk 


nidkUj + nidjfik - njUinmOmUk 


+ 4 e jik 


nid k rii + nidiUk - 


(17) 


The issue at hand is what a suitable replacement should 
be if we express the orientational order using a Q-tensor 
rather than the director field. The usual approach to ob¬ 
taining such a replacement is to adopt algebraic methods 
rather than geometric ones 1105. : one constructs a list of 
potential expressions and selects on the basis of how they 
reduce when the Q-tensor is written in terms of the di¬ 
rector field. Following this procedure, we find that the 
expression 


fib? — . ^ilk 


+ 4 e 3lk 


2Qi m dkQj m T QjmdkQl m 
+ Ql m djQk m Qjl^mQkn 
2 Ql m dkQ im T Qim dkQ i r 
+ Ql m diQk m QildmQkrr 


(18) 


reduces to n ^ when the Q-tensor is uniaxial, with con¬ 
stant magnitude, Qij oc nirij — \5ij- We subtract its 
trace and adopt it as an appropriate Q-tensor analogue 
of n. In simulations, the umbilic lines can then be iden¬ 
tified using isosurfaces where the norm of n drops below 
a threshold value. 

Illustrations are given in Fig. [ 6 ] for some standard 
blue phase textures. The umbilics coincide with struc¬ 
tural motifs that have been identified previously, either 
through symmetry considerations m or by taking a 
singular value decomposition of the Q-tensor [69]. These 
motifs are the axes of double twist cylinders [ 6 ] I. Two 
remarks are in order concerning this. First, the pattern 
of umbilic lines in BPII (see Fig. | 6 | does not correspond 
to the array of double twist cylinders traditionally de¬ 
picted in the literature. In the traditional depiction the 

























FIG. 6: Umbilic lines in the blue phases. The umbilics are 
indicated by blue isosurfaces, while the network of disclination 
lines is shown by milk chocolate isosurfaces. The textures BPI 
and BPII are the cubic structures observed experimentally. 
08+ is a cubic structure with the same space group as BPI 
but an interchange between disclination lines and umbilics. 
BPX is a tetragonal texture observed in an applied electric 
field. The arrangement of disclinations and umbilics is the 
same as in O s+ , although the symmetry is different. 

double twist cylinders are taken along the 2 -fold axes ly¬ 
ing in the faces of the conventional unit cell and, indeed, 
inspection by eye reveals these lines to have the same 
general structure as model double twist cylinders [6 (see 
Fig. 0 - However, in a Q-tensor description of BPII these 
are not degeneracies where the local structure is axisym- 
metric, as in the standard idealised picture of a double 
twist cylinder; rather they are 2 -fold lines where the ex¬ 
tensor is maximally biaxial piiwinns] . Second, it is 
perhaps worth emphasising that the umbilics identified in 
all other blue phase textures also differ from the standard 
idealised picture of a double twist cylinder. This is be¬ 
cause the idealised axisymmetry of the latter constrains 
it to be a degenerate winding +2 umbilic, whereas those 
that are found in periodic textures are generic umbilics, 
typically with a — 1/2 profile for the eigenvectors of n, 
i.e. a star umbilic. 

We have given here only cursory illustrations for peri¬ 
odic blue phases with crystalline space groups. It would 
be of interest to also study the structure of the amor¬ 
phous BPIII in terms of umbilics, as well as their be¬ 
haviour in transitions and rheology. The Q-tensor ver¬ 
sion of the linear transformation n that we introduced 
here was motivated on algebraic grounds rather than the 
geometrical considerations that the rest of our work is 
founded upon. It is clear that geometrical considerations 


can be developed for traceless symmetric tensors, paral¬ 
lel to our description of orientational order. It would be 
interesting to do so and compare the outcome with the 
algebraic construction we have given. 

In the remainder of this paper we return to the setting 
of vector orientational order. 

LOCAL PROFILES OF UMBILIC LINES 

Generic umbilic lines carry a canonical orientation 
through the vector T>, defined as 

B i = (n, cFA) = (A, A )6 ij Aj, (19) 

which assigns a tangent vector to the umbilic, provided 
it is generic, as follows. If, in a local trivialisation, A = 
siei ± s 2 e 2 then an umbilic is defined as the intersection 
of the two surfaces si = 0 and s 2 = 0 and so its tangent 
vector will be Vsi x Vs 2 . Note that while the values of s i 
and 52 depend on the choice of basis, the tangent vector 
does not. In this local form we can write 

B = s 2 — 52 V 51 T 2(<s^ T ^ 2)^5 (20) 

so that on an umbilic we have V x B = 2Vsi x V «2 
and hence the vector Vx5 points along umbilic lines, 
canonically orienting them. As indicated in Fig. [3] this 
orientation is compatible with that induced by the circu¬ 
lation in A previously discussed. As alluded to in Figs.[l] 
and [4j the discrepancy between the type of profile and 
index of an umbilic depends on the relative orientation 
between n and the umbilic line. To measure whether 
the eigenvectors of a generic umbilic have winding ± 1/2 
one should integrate A along a contour 7 chosen such 
that its orientation matches the orientation imparted to 
£ by n. One can then show that the local profile is ± 1/2 
depending on 

Sgn(n • V x T>), ( 21 ) 

the analysis being identical to that given by Berry for C 
lines in electromagnetic fields m- 

This quantity n • V x B can be used to further char¬ 
acterise the local structure of umbilic lines. Along an 
umbilic A vanishes and generically we can expect it to 
vanish linearly. Its local profile is then given by the form 
of the linear order terms in A in the immediate vicinity 
of the umbilic. Equivalently, generically we can expect 
the gradients of A to be non-zero on the umbilic and the 
local profile is then given by the structure of VA eval¬ 
uated on the umbilic. As usual, we split the gradients 
into the components parallel to n and those perpendic¬ 
ular to it, and consider only the latter. The non-zero 
part of these transverse gradients, evaluated on the um¬ 
bilics, takes values in ® E , so it is this part that we 
need to focus on. A section of this bundle characteris¬ 
ing the local profile of the umbilic can be constructed 
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from VA as follows. If Ie is the projection onto E in 
TM 3 * 0 TM 3 then we define the differential operator 
as V^A = I^VAI#. By construction V^A is a section of 
the bundle £* ®E. Under the action of 50(2) the bundle 
£*0T splits into two rank 2 subbundles, T + ®T _ , where 
T + is the spin 1 bundle corresponding to +1/2 profiles 
and F~ is the spin 3 bundle corresponding to —1/2 pro¬ 
files, and we decompose V^A according to this splitting 
as 


V^A = V^A+ + V*A“. (22) 

A short calculation shows that n • V x B = |V^A + | 2 — 
|V^A“| 2 so that if |V^A + | > |V^A - | the umbilic has a 
local +1/2 profile, while if the converse holds the profile 
is of type —1/2. 

This analysis allows us to give a description of the 
space of local profiles for a generic umbilic. Since V^A|^ 
is non-zero for a generic umbilic, but otherwise arbitrary, 
the space of local profiles has the homotopy type of the 
three-sphere. However, there is more structure to it than 
that. The equality |V^A + | = |V^A“| defines a Clifford 
torus separating the three-sphere into two solid tori cor¬ 
responding to the different profiles. A natural constraint 
on an umbilic loop is for its profile to be of constant type, 
i.e. for n • V x B to be of constant sign or the profile 
to stay within a single solid torus. We call an umbilic 
loop which satisfies this constraint transverse, as its tan¬ 
gent vector is always transverse to the planes of £. This 
is not an unphysical constraint, transverse umbilic lines 
are commonly found in various chiral structures - typi¬ 
cal A lines in cholesterics have a local structure in which 
n is either parallel or antiparallel with the tangent vec¬ 
tor to an umbilic. Figure [7] shows the transverse umbilic 
loop found in toron excitations in frustrated cholester¬ 
ics [32| l3i ] . 

The condition of an umbilic being transverse there¬ 
fore translates as one of the two profiles having a higher 
‘weight’ along the entire length of the umbilic. If this 
is the case then we have |V^A ± | 2 ^ 0 with the sign de¬ 
pending on whether the umbilic is positively or negatively 
transverse. Suppose, for concreteness, that the umbilic 
is positively transverse. Then, in a local trivialisation, 
we will have V^A + = t\ f+ + with t\ + 1\ / 0, and 
the space of positively transverse umbilics has the homo¬ 
topy type of a circle. The variation of the profile around 
a closed, transverse umbilic loop confers another integer 
winding number. This number depends on how the ba¬ 
sis vectors of the trivialisation of T + wind around the 
umbilic. The ambiguity can be removed by demanding 
that each of the basis vectors in the trivialisation does 
not link with the umbilic. As we will show in the next 
section, this winding number can be interpreted as twice 
a self-linking number for the umbilic loop, and is related 
to a relative Euler class of E. 



FIG. 7: Umbilic lines in a toron. (a) Simulation results show¬ 
ing point defects and characteristic umbilic loop in a toron. 
The umbilic loop is transverse and has zero self-linking num¬ 
ber. The point defects have opposite charge, due to lattice 
effects the simulation does not resolve the umbilic lines con¬ 
necting them, though they are clear in the topology of the 
eigenvector field, (b) Diagram illustrating the umbilic struc¬ 
ture of the toron in S 3 when the boundary of the cell is col¬ 
lapsed to a point, note that the umbilic loops connecting the 
point defects each have strength s = 2, for a total of 4, as 
required by ( fl3] ) . (c) Simulation results showing director field 
for the toromfd) Contour plot showing the value of |A|. (e) 
Eigenvector field of n. Note that the umbilic loop is generic 
and has a monstar profile. 


UMBILIC LOOPS 


The umbilics of three-dimensional orientationally or¬ 
dered materials are extended line-like objects, however 
most of the properties we have described so far can be 
associated to isolated points along these one-dimensional 
curves. For example the points of transverse intersec¬ 
tion with an appropriate surface (or slice through the 
material) count numbers of Skyrmions, merons, or A de¬ 
fects. It is clear, however, that these locally measurable 
objects do not represent all the information encoded by 
umbilic lines; if one has a closed umbilic loop then, as dis¬ 
cussed above, the local proprties at each point must stitch 
together consistently along the loop conveying both ge¬ 
ometric and topological information about the orienta¬ 
tional order. 

Umbilic loops have been observed in experimental sys¬ 
tems mmi Em, and a motivating set of examples are the 
so-called torons [32 , illustrated in Fig. [7J These soliton 
structures can be generated in a frustrated cholesteric 
cell, with the director held vertical along the top and 
bottom. As shown in Fig. [7|[a), they display a charac- 
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terisitc umbilic loop, in the form of a circle sitting in the 
midplane of the cell. This loop is of generic type, with 
the eigenvectors of II having a monstar profile. 

As with local properties of umbilic lines, much of the 
structure of loops is encoded in the 1-form A. Previ¬ 
ously we have obtained local information about umbilics 
by integrating A around a meridian. The first step to un¬ 
derstanding their global structure is to integrate A along 
a longitude. Since A is singular on an umbilic, we must 
define this integral as a limit. If Ui is an umbilic, then let 
U[ be a curve close to Ui that has zero linking number 
with it. The integral of A over U[ is well defined and 
if all the umbilics in the system are closed loops we can 
take the limit U[ Ui to obtain 



2ir^s j Lk(U i ,U j ) + j Q e , 


(23) 


where E is an orientable surface bounded by Ui. The 
requirement that U[ have zero linking with Ui ensures 
that the limit is well defined but in general the re¬ 
sult depends upon this choice of framing. The quantity 
JA SjLk(Ui, Uj) = e(Ui) in (23) is an integer associated 
to the umbilic line, the relative Euler class, and counts 
the number of umbilics linking Ui. There are two natural 
questions to ask of this quantity. The first is whether it 
corresponds to any local geometric structure of the um¬ 
bilic line; the second is whether it relates to a global 
property of n. On both counts the answer is yes. Geo¬ 
metrically (23) is related to both a self-linking number 
and the winding number defined in the previous section. 
Globally we will show that it is related to the Hopf in¬ 
variant of n. 


To see this suppose U is a transverse umbilic loop, then 
let U be a loop close to U of zero linking number with 
U. Evaluate one of the eigenvectors of A, say, along 
U and then push U along this eigenvector field to form 
a new loop, U' . The linking number Lk (U,U f ) is then 
taken to define the self-linking number of the transverse 
umbilic, SI (U). This construction is shown in Fig. [8j 
Note that since pj is a line field, SI (U) may be a half¬ 
integer corresponding to p being non-orientable along the 
longitude of U. 

With this notion in hand the correspondence between 
the self-linking number and the winding number coming 
from the variation of the local profile defined in the pre¬ 
vious section can be established. Choose a trivialisation 
{d!,d 2 } of £ on a neighbourhood of U such that the link¬ 
ing number of U pushed off along di (say) is zero then 
the winding of V^A + with respect to the chosen triv¬ 
ialisation will give twice the self-linking number. Note 
that the geometric definition of the self-linking number 
did not depend on the zero of A being of linear order. 
The relationship between between the self-linking num¬ 
ber and the relative Euler class of a transverse umbilic is 



FIG. 8: The self-linking number for transverse umbilic loops. 
The red curve U' is formed by pushing U along an eigenvector 
field pa of A. The self-linking number is then computed 
as the linking number of U with U '. The umbilic must be 
transverse to ensure that U' does not cross U . Because the 
eigenvectors are degenerate on the umbilic, one evaluates pa 
on a line close to U , U, that has zero self-linking number 
with U. Such a line can be found by drawing an orientable 
surface whose boundary is U and pushing U into that surface. 
Note that because eigenvectors are line fields, the self-linking 
number can be a half-integer, corresponding to p being non- 
orientable along the longitude of U. 


then |110| 

S1(C) = SI (CO + \e(U). (24) 


There is a new term in this equation, SI (I/), the trans¬ 
verse self-linking number of U. SI (U) is a geometric quan¬ 
tity, preserved under deformations that keep U trans¬ 
verse. A similar notion is studied in the context of con¬ 
tact topology m where it provides useful geometric 
information, but we do not know of a general discussion 
in the case of arbitrary n. 

We now turn to the relationship between the linking 
of umbilics and global properties of n. Associated to the 
1-form A is the Chern-Simons 3-form, A A cL4, and it is 
natural to ask whether the integral over the complement 
of W, 


/ 


A A <A4, 


(25) 


gives some information about umbilic lines and n. This 
integral can be compared to the similar helicity integral 
for a fluid flow and the Abelian Chern-Simons action. 
Both these quantities are related to linking numbers, of 
vortex lines in the case of helicity [64lf66] and Wilson 
loops in Chern-Simons theory (112| II 13] and, as we will 
show, it is no different for umbilics, (25) connects the 
linking numbers of the umbilic loops with the Hopf in¬ 
variant of n. To proceed we will assume that the material 
domain is S' 3 , with no defects and that the zero set of A 
is one-dimensional. By a compactification, this assump¬ 
tion that the domain is S 3 is equivalent to examining a 
system in M 3 for which lirm^oo n(x) = no, a constant. 
An example of such a system is illustrated by the Hopf 
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FIG. 9: Umbilic lines in the Hopf fibration. Top : Umbilic 
lines given as isosurfaces of | A| for a simulated Hopf fibration. 
The two lines in the mid plane of the cell are both generic, 
with the same orientation. For symmetry reasons the central 
umbilic is not generic and has s = 2. On the boundary of 
the cell n is uniform, which means that the Hopf invariant is 
well-defined. Note that the system is totally umbilic on the 
boundary. If you compactify the domain by collapsing the 
boundary to a point, then one can interpret the central ver¬ 
tical umbilic line as passing through this point. Bottom Left : 
director field for the Hopf texture. Bottom Right : Schematic 
of the relationship between umbilic lines in S' 3 . Two generic 
s = 1 umbilic lines each link with one s = 2 umbilic line that 
passes through the point at infinity. 


fibration, which has been experimentally realised in frus¬ 
trated cholesterics [34] . illustrated in Fig.[9j and the afor- 
mentioned toron. In both these systems n is constant on 
the boundary of the cell. Because of this A = 0 on the 
boundary and one should consider these systems as hav¬ 
ing an umbilic line passing through the point ‘at infinity’, 
which in Figs. [7] and [9] is just the boundary of the system. 

Topologically the configuration in Fig. [9] is unusual, 
it is non-singular but not homotopically equivalent to a 
constant. The difference is that it has a non-zero Hopf 
invariant. This integer invariant distinguishes topologi¬ 
cally distinct non-singular textures in M 3 which are con¬ 
strained to tend to a constant at infinity and is typically 
visualised in terms of linking of pre-images of two orien¬ 
tations using the Pontryagin-Thom construction |H]. In 
our setting however one can give a slight re-interpretation 
of the famous Whitehead integral formula m, similar 
to that given by Arnol’d |66 to show that the Hopf in¬ 
variant is also the Chern-Simons invariant of the bundle 
£. In our interpretation, if one takes oo to be a connec¬ 
tion form for £, then duo is the curvature of £ and the 


(47r) : 


/■ 

Js 3 


oo A duo 


(26) 


computes the Hopf invariant of n. 

Now we turn back to our integral ( [25| . As discussed, 
we assume the material domain is S 3 (or M 3 with uniform 
boundary conditions). This means that we can give a 
global trivialisation of £ with connection oo and our 1- 
form A can be written as A = d6 + 2oo. This allows us to 
write 

/ AAgL4 = 4 / oo A duo+ 2 dOAdoo (27) 

Js 3 -u Js 3 -u Js 3 -u 

or, since U has measure zero 

I A/\dA = 4(47r) 2 iJ + 2 [ dOAdoo (28) 

Js 3 -u Js 3 -u 

where H is the Hopf invariant of n. Now we need to 
evaluate the integral of dO A duo over S 3 —U. We will first 
do a simpler integral, over S 3 — N(U), where N(U) is a 
neighbourhood of the umbilic lines. By Stokes’ theorem 
this can be reduced to an integral over the boundary of 
N(U) as 


L 


dO A duo = 


S 3 —N(U) 


?/. 


lo A d6. 


dN(Ui) 


(29) 


The boundary of the neighbourhood of each umbilic line 
is a torus and because our domain is S 3 we can choose 
a preferred longitude for each torus that has zero linking 
with the umbilic at its centre. This allows us to specify 
a decomposition of each boundary component into S]^ x 
S\. If we choose N(U) small enough, then along each 
meridian of the torus oo will be approximately constant 
and equal to its value on the umbilic line. The integral 
over each boundary component can then be written as: 


?/. 


oo A dO = 


dN(Ui ) 


;/ -*(/ 

J si \Jsi 


de* 


(30) 


where we have pulled the forms back to each S l . 

The integral of dO is trivial, it evaluates to 5^277, where 
Si is the strength of the i th umbilic. Finally, we take a 
limit as N(U) —» U and obtain a sum of Wilson loop 
integrals 


lim 

N(U)^U A 


/ oo A dO = 2tt Si / 

„• JdN(Ui) , JUi 


00. 


(31) 


Now from the Gauss-Bonnet-Chern theorem (13) we can 


relate inter gals of oo along umbilics to integrals of A. 
Putting this together with the above equations we ob¬ 
tain 


A 


- 2 — [ A A <L4 - T W s . [ 

(4?r) 2 Js 3 _ u 8tt A 

= 4F-N S , Sj L k (f/ ! ,^) (32) 
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This equation relates analytic data, the integral of A Ad A 
and Wilson loop integrals, to topological data, namely 
the linking number of the umbilic lines and the Hopf in¬ 
variant. We hope that this relationship will allow for 
an interpretation of the creation and destruction of Hopf 
solitons in terms of umbilic lines and lead to connections 
between the dynamics of umbilic lines and of vortex lines 
in fluids. It is also tempting to speculate on the use um- 
bilic loops as generators in a Floer theory of orientational 
order. 
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